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Abstract — Optimal power flow (OPF) is considered for micro- 
grids, with the objective of minimizing either the power distribu- 
tion losses, or, the cost of power drawn from the substation and 
supplied by distributed generation (DG) units, while effecting 
voltage regulation. The microgrid is unbalanced, due to unequal 
loads in each phase and non-equilateral conductor spacings on 
the distribution lines. Similar to OPF formulations for balanced 
systems, the considered OPF problem is nonconvex. Nevertheless, 
a semideflnite programming (SDP) relaxation technique is advo- 
cated to obtain a convex problem solvable in polynomial-time 
complexity. Enticingly, numerical tests demonstrate the ability of 
the proposed method to attain the globally optimal solution of the 
original nonconvex OPF. To ensure scalability with respect to the 
number of nodes, robustness to isolated communication outages, 
and data privacy and integrity, the proposed SDP is solved in 
a distributed fashion by resorting to the alternating direction 
method of multipliers. The resulting algorithm entails iterative 
message-passing among groups of consumers and guarantees 
faster convergence compared to competing alternatives. 

Index Terms — Microgrids, distribution feeders, optimal power 
flow, semideflnite relaxation, distributed computation, alternating 
direction method of multipliers. 



I. Introduction 

Microgrids are portions of a power distribution network 
located downstream of the distribution substation that supply 
a number of industrial and residential loads, and may include 
distributed generation (DG) and energy storage devices [1]. 
A microgrid can operate in either grid-connected, islanded, 
or hybrid modes. Deployment of microgrids promises drastic 
performance enhancement of the distribution grid in terms of 
efficiency and stability, along with increased network scalabil- 
ity and resilience to outages. 

Besides bringing power generation closer to the end user, 
DG units offer environment-friendly advantages over conven- 
tional generation [1], may provide ancillary services such 
as reactive and harmonic compensation [2], [3], and enable 
DG owners to actively participate in grid operations through 
supply contracts and pricing schemes. On the other hand, their 
operation must be carefully controlled in order to prevent 
abrupt voltage fluctuations, which stem from the well known 
sensitivity of voltages to variations of power injections, node 
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over- and under- voltages [4], and drops of the power factor 
at the distribution substation [5]. To this end, optimal power 
flow (OPF) approaches are increasingly advocated also in this 
context to ensure efficient operation of smart microgrids and 
effect strict voltage regulation [3], [5], [6], [7]. 

OPF problems are deemed challenging because they re- 
quire solving nonconvex problems. Nonconvexity stems from 
the nonlinear relationship between voltages and the complex 
powers demanded or injected at the nodes. In the context 
of transmission networks, the Newton-Raphson method [8] 
is traditionally employed to obtain a possibly suboptimal 
solution of these nonconvex problems. In distribution systems 
however, its convergence is challenged by the high resistance- 
to-reactance ratio of distribution lines. Alternative approaches 
include sequential quadratic optimization, steepest descent- 
based methods [5], fuzzy dynamic programming [9], and 
particle swarm optimization [7]. However, these methods 
generally return suboptimal load flow solutions, and may be 
computationally cumbersome [10]. To alleviate these concerns, 
a relaxed semideflnite programming (SDP) reformulation of 
the OPF problem for balanced transmission systems was 
proposed in [11] and [12], where global optimality can be 
assessed by checking the rank of the obtained voltage-related 
matrix. The relaxed SDP approach was extended to balanced 
distribution systems in [13] and [6]. Notably, for networks with 
a tree topology, [6] and [13] established sufficient conditions 
under which a globally optimal solution is attainable provided 
the original OPF problem is feasible. 

Distribution networks are inherently unbalanced because: 
i) unequal single-phase loads must be served, and ii) non- 
equilateral conductor spacings of three-phase line segments 
are involved [14]. Further, single-phase DG units may worsen 
the network imbalance. As a consequence, optimization ap- 
proaches can not rely on single-phase equivalent models as in 
e.g., [3], [5], [6], [7], [13]. For the unbalanced setup, an OPF 
framework was proposed in [15], where commercial solvers 
of nonlinear programs were adopted, and in [16], where quasi - 
Newton methods were utilized in conjunction with load flow 
solvers. However, since these methods are inherently related to 
gradient descent solvers of nonconvex programs, they inherit 
the limitations of being sensitive to initialization, and do not 
guarantee global optimality of their solutions. 

The first contribution of the present paper consists in perme- 
ating the benefits of SDP relaxation techniques [17] to OPF 
problems for microgrids operating in an unbalanced setup. 
This optimization tool has three main advantages: ;) it offers 
the potential of finding the globally optimal solution; ii) its 
worst-case computational complexity is quantifiable; and, Hi) 
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it can accommodate additional thermal and quality-of-power 
constraints without exacerbating the problem complexity [17]. 
Global optimality not only reduces power distribution losses, 
but also leads to higher monetary savings compared to subop- 
timal OPF solutions. 

The OPF problem is solved by a microgrid energy manager 
(MEM), which cooperates with local controllers (LCs) located 
throughout the network. Microgrids can vary in scope, size, 
and ownership [18]. For those of medium and large size (of 
a distribution feeder), solving the OPF problem centrally at 
the MEM may become computationally prohibitive [10]. In 
fact, interior point SDP solvers do not generally scale well 
with the problem size [6], [19], [17]. For a real-time network 
management, it is generally required to find a new network 
operational setup rapidly (e.g., in a few seconds or minutes) 
in order to promptly respond to abrupt load variations and to 
cope with the intermittent power generation that is typical of 
renewable-based DG units. It is then of paramount importance 
to solve the SDP-based OPF problem in a distributed manner, 
by decomposing the main problem into multiple sub-instances 
that can be solved efficiently and in parallel. A distributed 
algorithm is desirable also when (a group of) customers do 
not share data with the MEM due to privacy concerns, or 
because they wish to manage autonomously their DG units in 
order to pursue individual economic interests [18]. Finally, 
a distributed algorithms involves a modest communication 
overhead compared to its centralized counterpart, as it does 
not require to pool line, generator, and load data at the MEM, 
and subsequently disseminate the OPF solution. 

Decentralized OPF approaches were first proposed in [20], 
[21], where multi-utility transmission systems were partitioned 
autonomously managed areas. Augmented Lagrangian meth- 
ods were employed to decompose the overall OPF problem 
in per-area instances. A similar approach was followed by 
e.g., [22], [23] (see also references therein), where standard 
Lagrangian approaches were utilized in conjunction with 
Newton methods. Solving SDP in a distributed fashion is 
challenging due to the couplings of local voltage-related 
matrices enforced by the positive semidefinite constraint of 
the global voltage matrix. Results related to positive semidef- 
inite matrix completion [24] were leveraged in [6], [19] to 
develop a distributed OPF algorithms for balanced networks 
via dual decomposition. Using the results of [24], and tapping 
into the powerful alternating direction method of multipliers 
(ADMM) [25, Sec. 3.4], a distributed optimization problem 
for unbalanced microgrids is formulated here, where each LC 
solves an optimization sub-problem, and then exchanges sim- 
ple messages with its neighboring LCs. Compared to the dual 
decomposition schemes of [6], [19], the proposed ADMM- 
based approach offers a markedly improved convergence. 

In the OPF context, augmented Lagrangian methods (re- 
lated to ADMM) were first used in [20], [21] to develop a 
decentralized optimization scheme for (balanced) transmission 
networks. Off-the-shelf schemes were used to solve the OPF 
sub-problem associated with each sub-network. More recently, 
ADMM was advocated in [26] for distributed state estimation, 
and in [27] for distributed multi-period OPF in balanced sys- 
tems. Here, the approaches of [20], [21], [27] are considerably 



broadened by considering unbalanced distribution networks, 
and by employing the ADMM to devise a distributed SDP 
solver. 

The rest of the paper is organized as follows. Section [II] 
recapitulates the OPF formulation for microgrids, and Sec- 
tion [III] develops its centralized SDP solver. The distributed 
algorithm is presented in Section [IV] while numerical tests 
are reported in Section [V] Finally, concluding remarks are 
provided in Section 

II. Modeling and problem formulation 

Consider a microgrid comprising N nodes collected in the 
set Af := {0, 1, . . . , N}, and overhead or underground lines 
represented by the set of edges^ £ := {(m,n)} C Af x Af. 
Assume that the network has a tree topology, and let node 
represent the point of common coupling (PCC), taken to be 
the distribution substation. Define as V mn Q {a>mn, b m n, c mn } 
and V n C {a n , b n , c n } the phases of line (m, n) G £ and node 
n E Af, respectively. Let £ C be the complex line-to- 
ground voltage at node n S Af of phase <j> 6 P„, and 1% € C 
the current injected at the same node and phase. As usual, 
the voltages vo := [Vq, V h , Vq] 7 ' are taken as reference for 
the phasorial representation. Lines (m, n) <E £ are modeled 
as 7r-equivalent components [14, Ch. 6], and the \V m n\ x 
\P m n\ phase impedance and shunt admittance matrices are 
denoted as Z mn e C^^^ and Y$ n e C^™ l x l^ I, 
respectively. Three- or single-phase transformers (if any) are 
modeled as series components with transmission parameters 
that depend on the connection type [14, Ch. 8], [15]. 

Per phase e V n , let Pf in and Q^ L denote the active 
and reactive powers demanded by a wye-connected load at 
the bus n. Capacitor banks are usually mounted at some 
nodes to provide reactive power support. Let denote the 
susceptance of a capacitor connected at node n and phase <j>. 
Finally, suppose that S DG units are located at a subset of 
nodes S C Af, and let Pq and denote the active and 

reactive powers supplied by unit s E S. For conventional DG 
units, such as diesel generators and fuel cells, the supplied 
powers can be controlled, and they will be variables of the 
OPF problem; on the other hand, Pq and Qq represent 
committed powers for renewable-based sources. 

Given the demanded loads {P£ n , n }, the goal is to 
select a feasible set of voltages {V^f}, currents and 
powers supplied by conventional DG units {Pq s ,Qq s } so 
that the steady-state operation of the microgrid is optimal in 
a well defined sense. To this end, one of the following two 
objectives is usually pursued. 

1 Notation: Upper (lower) boldface letters will be used for matrices (column 
vectors); (•) 7 " for transposition; (■)* complex-conjugate; and, complex- 
conjugate transposition; 3?{-} denotes the real part, and the imaginary 
part; j := the imaginary unit. Tr(-) the matrix trace; rank(-) the matrix 
rank; | • | denotes the magnitude of a number or the cardinality of a set; and 
|| ■ \\p stands for Frobenius norm. Given a vector v and a matrix V, [v]-p 
denotes a \V\ X 1 sub- vector containing entries of v indexed by the set V, and 
[V]-pj p, the |"Pi| X \Py\ sub-matrix with row and column indexes described 
by V\ and "Pi. Finally, OmxN and liVfxiV denotes M X TV matrices with 
all zeroes and ones, respectively. 

2 The symbols defined throughout the paper are recapitulated in Table [J 
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i) Minimization of power losses. The active power flowing 
on line (m,n) G £ is P m ^ n ■= E^, e -p ran ^(^Ln)*' where 
Im^t-n i s me complex current from node m to node n on phase 
G T^mn- Thus, the cost to be minimized here is [6] 



d(V) 



E 

(m,n)G£ 



(i) 



where V := },{#}, {P£ n ,Q£, n }} collects a11 the 

steady-state variables. 

ii) Minimization of the cost of supplied power. Letting the cost 
of power drawn from the PCC be denoted by Co > 0, and the 
one incurred by the use of DG unit s G S as c s > 0, one can 
minimize the cost of supplied power [5] 

E c 

ses <pev s 

Notice that ((TJ and (f2]i are equivalent when Co = 1 and c s = 1 
for all s G S. 

Based on (fTJ-ffJJi, the following OPF problem is considered: 



C 2 (V) 



(2) 



(PI) 
s.t. Vf(lf) 



min C m (V) 
v 



(3a) 



P 



it = 



V 4>£V s ,seS (3b) 

-^.n-iOin + i^nlKfl 3 . 

V^P„,ne Af\S (3c) 
£ [(iyw +z- 1 n )[v n ]^ ron 



T/IJ 



_ Z mn [v m ]p„ 

< iKf I < vr\ 



{'!>} 



pg™ < fg, a < PS™, 



< Q%, s < Q 



max 



V <f> G P n , n G M (3d) 

V G P„, n&Af (3e) 

V 4> G P s , s G S (3f) 

V cj) £P S , s G <S (3g) 



where V^" m and V^ nax in (I3eb are given minimum and 
maximum utilization and service voltages; J\f n := {j\(n,j) G 
£}; (T3fb — (|3g]> are box constraints for the power supplied by 
conventional DG units, and m G {1,2} depending on the 
chosen cost. 

Similar to OPF variants for transmission networks and 
balanced distribution networks, (PI) is a nonlinear nonconvex 
problem because of the load flow equations (|3bb-(l3cT> as well 
as the constraints (|3e] >. In the next section, an equivalent 
reformulation of (PI) will be introduced, and its solution will 
be pursued using the SDP relaxation technique. 

III. SDP-Based Centralized Solution 
Consider the ~^2 n& j^\P n \ x 1 complex vectors v 



f > v i j ■ 



,v^] r and i 



;7 :T 



• ) 1 N 



ijj] 7 ', with v„ and 



i„ the \P n \ x 1 complex vectors collecting voltages {V^} c j >( z-p n 
and currents {In}<bev n per node n G N. Voltages and injected 
currents abide by Ohm's law i = Yv, where Y is a symmetric 
block matrix of dimensions J2neM \Pn \ x J2neAf |P«I' wnose 
entries are given by: 

i) matrix — occupying the \P mn \ x \V mn \ off-diagonal 
block corresponding to line (m, n) G £; and, 



TABLE I 
Nomenclature 



£ 
S 

M n 

'Pn 
v n 

jo 

1 n 

p<t> 

L,n 
p<l> 
^G,n 

4> 

fc,n 

Pm — > n 

v < 3 ) 

1 mn 

ymin ymax 
pmin ( /">min \ 
pmax /nmax \ 



«&.») 



AW 



Set collecting the nodes of the microgrid 

Set collecting the distribution lines of the microgrid 

Subset of nodes featuring DG 

Set of neighboring nodes of n S M 

Set of phases at node n 

Set of phases of line (m, n) 

Complex line-to-ground voltage at phase <j> of node n 
Complex current injected at phase <p of node n 
Complex current on phase (f> of line (m, n) 
Active (reactive) load demanded at node n on phase <f> 
Active (reactive) power supplied at node n on phase cj> 
Susceptance of a capacitor at node n and phase <f> 
Active power exiting node m on line (m, n) 
Phase impedance matrix of line (m, n) 
Shunt admittance matrix of line (m, n) 
Utilization and service voltage magnitude limits 
Minimum active (reactive) power supplied at node n 
Maximum active (reactive) power supplied at node n 
Set of nodes forming the area I of the microgrid 
Extended area I comprising the nodes in A^' and the 
nodes of different areas connected to A^ ' by a line 
Set of neighboring areas of AW 



ii) the \P n \ x \P n \ diagonal block corresponding to node 
ii£jV with U m := {n\(m, n) G £ } 



v n ,v„ 



E 



Yl s l + Y„ 



7-1 
mn 
-1 



if P-n = 



Zj n and [Y, n „] 



mn\V n \V nm ,V n \P n 



(4) 



otherwise 
= 



where Y mn 

[Ymn]'Pnm.,'Pnm 

(Ymn is formed likewise). 

The next step consists in expressing the active and reactive 
powers injected per node, active powers flowing on the lines, 
and voltage magnitudes, as linear functions of the outer- 
product matrix V := vv M . To this end, define the following 
admittance-related matrix per node n and phase (f> 

Yt := et(et) r Y 



where := 



l°Tvo\> 



' u |p„-il 



b.T 



\> n+l\ 



;0T 



(5) 
IT 



and {ef, j^ep^ denotes the canonical basis of Rl 73 "!, 
and let the \P mn \ x Y^neM^^ matrices A ™^« 

respectively. Then, a linear model in V can be established 
using the following lemma (see also [28] and [12]). 



and B,„ 

V 1 n 

B„, := 



be defined as A m _> 

_ry-l 

«n|x£r=m+l " r 

[0, 



Lemma 1: For the Hermitian matrices 



P,n ■- 



3 



:= J -(Y£-(Yt)*) 



V,n 



eS(e*) 



T 



H 



B, 



R w A 



(6a) 

(6b) 
(6c) 
(6d) 
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voltage magnitudes and active as well as reactive powers are 
linearly related with V as 



\v*\ 2 

p<t> _ -p<t> 

r G,n r L,n 

Tr(#£ )n V) = Qt. n Qt n + ^,„Tr(*t,„V) 
Tr(* m _v„V) = P m ^ n 



■ V,n 

Tr(*t„V) 



with P* tn = 



(7a) 
(7b) 
(7c) 
(7d) 

= if 

□ 



S,„ = for n e Af\S, and y* „ 
capacitor banks are not present at node n. 
Proof. See the Appendix. 

Using Lemma Q] problem (PI) is equivalently reformulated 
as follows: 

(P2) mm (7 m (V) 

s.t. Tr(*£ n V) + P£ t „ = 0, V 0, V n e M\S (8a) 

Tr(** in V) + Q(„ - y* „Tr(** n V) = 0, 

V 0, Vn eA/AS (8b) 
Pg^ < Tr(** a V) + Pt s < P£™, V <j>, Vn G 5 (8c) 
< Tr(*^ >s V) + g*, s < QS™, V 0, Vn e 5 (8d) 
(K min )2 < Tr($t,„V) < (l/r x ) 2 , V 0, Vn e A" (8e) 
0, [VK,p =v v^ (8f) 
rank(V) = 1 (8g) 

where the costs C\ (V) and Ci (V) are re-expressed as 

^i( V ) = Yu (Tr(* m ->„V)+Tr(* n _> m V)) (9a) 

(m,n)££ 

Problem (P2) is still nonconvex because of the rank-1 



constraint (8gi. Nevertheless, the SDP relaxation technique, 
which amounts to dropping the rank constraint [17], can be 
leveraged to obtain the following convex problem: 



(P3) 



min C m (V) 
v 



(10) 



.t. dHD — (ED . 



Clearly, if the optimal solution V opt of (P3) has rank 1, then it 
is a globally optimal solution also for the nonconvex problem 
(P2). Further, since (PI) and (P2) are equivalent, there exists 



v opt v^ t , and the optimal costs 



a vector v opt so that V opt 
of (PI) and (P2) coincide at the optimum. This is formally 
summarized next. 

Proposition 1: Let V opt denote the optimal solution of the 
SDP (P3), and assume that rank(V opt ) = 1. Then, the vector of 
line-to-ground voltages v opt := y/Xim, where Ai G R + is the 
unique non-zero eigenvalue of V opt and ui the corresponding 
eigenvector, is a globally optimal solution of (PI). □ 

The upshot of the proposed formulation is that a globally 
optimal solution of (P2) (and hence (PI)) can be obtained 
via standard interior-point solvers in polynomial time. In 
fact, the worst-case complexity of (P3) is on the order 



0( maX {iV c , £„ eA/ - \V n \] VX^v \ V n\ log(l/e)) for gen- 
eral purpose SDP solvers, with N c denoting the total number 
of constraints and e > a given solution accuracy [17]. 
Notice however that the sparsity of {3>p n , <I?q n , &y n } and 
the so-called chordal structure of a matrix can be exploited 
to obtain substantial computational savings; see e.g., [29]. In 
contrast, gradient descent-based solvers for nonconvex pro- 
grams, sequential quadratic programming, and particle swarm 
optimization, do not guarantee global optimality of the ob- 
tained solutions and are sensitive to initialization. Here, global 
optimality translates to lower distribution losses and increased 
monetary savings compared to sub-optimal OPF solutions. 

Since (P3) is a relaxed version of (P2), V opt could have rank 
greater than 1. In this case, rank reduction techniques may 
be employed to find a feasible rank-1 approximation of V opt 
provided it exists. For instance, the randomization technique 
offers a viable way to obtain a rank-1 approximation with 
quantifiable approximation error; see e.g., [17] and references 
therein. Albeit feasible for (P2), the resultant solution is 
generally suboptimum [17]. For balanced tree distribution net- 
works, [13] and [6] established conditions under which a rank- 
1 solution is attainable provided the original OPF problem 
is feasible. Unfortunately, when the tree power network is 
unbalanced, the results of [13] and [6] no longer apply, as 
explained in the ensuing Section IIII-BI However, an intuitive 
argument will be provided in Section IIII-BI to explain why a 
rank-1 solution is expected even in the unbalanced setup. But 
first, SDP-consistent constraints on line flows are derived in 
Section IIII-AI and a remark is provided. 

Remark 1. Step-down or in-line three- or single-phase trans- 
former banks (if any) can be accommodated in the formulated 
optimization problems by using their series component mod- 
els [14, Ch. 8], [15]. If a delta connection is employed on one 
side of the transformer, a small "dummy" resistance should 
be included between the primary and the secondary sides (one 
per phase) in order to ensure that the matrix V opt obtained by 
solving (P3) has rank 1; see also [12]. □ 

A. Constraints on line flows 

Constraints on the power dispelled on the distribution lines, 
or, on the line current magnitudes are generally adopted to 
protect conductors from overheating (which may eventually 
trigger an outage event). Using Lemma Q] it turns out that the 
real power dissipated on a line (m, n) £ £ can be limited by 
simply adding the constraint Tr(<fr m _>. n V) + Tr(3> n _>. m V) < 
AP m „ in (P3), for a given maximum power loss AP mn . 

Consider now the constraint < 7™^, with 

a given upper bound on the magnitude of lf 7ln . Aiming 
to an SDP-consistent reformulation of this constraint, let 
imn := [{ltn\V denote the \V mn \ x 1 vector collecting 
the complex currents flowing through line (m,n) E £, 
and notice that i mn is related to voltages v m and v„ as 
i mn = Z^ n ([v m ]-p mn - [v„]p ra J. Next, define the \V mn \ x 
S„ e .yv complex matrix 
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where Z™ n is a |P mn | x \P m \ matrix with elements 
[Z™J ?m „, Pm „ = Z-i and [iZn\v mn ,v m \v mn = 0; like- 
wise, Z" m has dimensions \P m n\ x |Pn|> and its entries are 
fe„]p„„,P nn = -Z^„ and [Z^ n ]p m „ >73 „\p m „ = 0. Thus, 
an SDP-compliant re-formulation of the constraint on the 
current magnitude is possible as follows. 

Lemma 2: Consider the Hermitian matrix 

<f>^ B n (e^ ) r B (12) 

where {e^j n }^ e p mn denotes the canonical basis of Rl Pm "L 
Then, constraint \I^ n \ < can be equivalently re- 

expressed as 

Tr{*L™V}<(/r) 2 - (B) 

□ 

Following similar steps, and using (l3dl i. constraints on the 
magnitude of the injected currents {1%} can be derived too. 

In unbalanced microgrids, it is of prime interest to protect 
from overheating also the neutral cable(s) of the distribution 
lines [14]. Towards this end, let Vmn denote the set of 
grounded neutral cables that are present on the line (to, n) E 
£. Further, let T mn the |Pron| x \V mn \ be the neutral transfor- 
mation matrix, which is obtained from the primitive impedance 
matrix of line (to, n) via Kron reduction [14, Sec. 4.1]. Thus, 
the neutral currents := [Imn, ■ • ■ , Imn n \ r are linearly 
related to the line currents i mn as imn = T m „i m „. It readily 
follows from Lemma [2] that the magnitude of the current on 
the neutral cables can be constrained in (P3) as 

Tr{*ff„„V} < (4^ max ) 2 , Vv>eP<$ (14) 

with — R w T H p^(p^h r T R lp (v Hthe 

Willi ^i mn i ■ — a mn ± mri K ™" 1 V^mn) J-ran D raiii i e mnj L11C 

canonical basis of M} v ^, and Imn' mdx a cap on the magnitude 
of 

B. The rank conundrum 

Sufficient conditions under which a rank-1 solution is al- 
ways obtained provided the original OPF problem is fea- 
sible were established in [13] and [6] for balanced distri- 
bution networks with a tree topology. Balanceness implies 
that equal single-phase loads are served, "L m l n = (g mn — 
jb mn )l\v nm \ and Y$ n = jb mn I\ Vnm \ for each line (m,n) <E 
£, where g m n,b mrl ,b mn > 0. To recapitulate the broad 
outline of the proofs in [13] and [6], assume for simplic- 
ity that the shunt admittances are all zero. Then, the total 
power flowing from node m to n is given by P m _>.„ = 
3|V^| 2 + 3\V*\\V*\(b mn sm(0t n ) - . 9m „cos«J), with 
®mn := @m ~ @t tne angle difference between voltages 
and V$ . Since the network is balanced, 0f nn is the 
same on each phase <f> £ Vmn- Fixing the voltage mag- 
nitudes, the region of the feasible powers (P m _>.„, P„_> TO ), 
which is denoted as J- mn , becomes an affine transforma- 
tion of the unit circles. Then, if one minimizes a strictly 
increasing function of the powers (P m _>. n , P n ->ro)> it follows 
that the Pareto front of J- mn and the one of the convex 
hull of T mn coincide if - tan" 1 (b nm /g nm ) < 6% m < 
tan _1 (6„ m /(7„ m ) [13]. Based on this observation, proving 




(a) Nonconvex problem. (b) Relaxed problem. 

Fig. 1. Line flow region and its Pareto optimal points. 




(a) Power injected at 3 nodes. (b) Injection at 2 and 3 (Pi = —5). 
Fig. 2. Feasible power injection region and its Pareto optimal points. 

that the Pareto regions of the feasible power injections at 
the nodes of the nonconvex OPF (the balanced counterpart 
of (P2)) and the relaxed SDP (the balanced counterpart of 
(P3)) amounts to showing that: i) the region of feasible 

powers T := {(P m -> )}(m,n)££\(P m ^ m , Prn^rn) € 

Fmn V (m, n) e £} is the Cartesian product of the regions 
{•7 r mn}(m,n)e£; an( l> U) the region of the injected powers is 
given by an affine transformation of T . Specifically, the first 
property i) follows from the fact that power flows on different 
lines are decoupled; that is, it is possible to modify the angle 
0f fln of a line (to, n), while preserving the angle difference 
Ott of any other line (fc, I) ^ (to, n). 

Suppose now that the off-diagonal elements of Z m „ are 
not zero; that is, ^ (g mn - jb mn )l\ Vnm \. The to- 

tal power flowing from node to to node n is given by 
P m ^n = 5i{v^Z~^(v m - v n )}, and it is now a function 
of {Q? nn }<f>£V mn > as we H as °f me angle differences {0 m ^ := 
&t-6 p m }<t>,pev m and {9*>p := 0*-0£}*,, 6 7v Different from 
the balanced case, the power flows {P m _>. n } are no longer 
decoupled across lines. In fact, it is impossible to adjust the 
angles {9 mn , Of^f , Of^ 9 } to obtain a new flow on line (to, n), 
without affecting the angle differences {9f , 9t ' p , Ot-f} and 
{0* k ,9fcP,e£' p } for one of the other lines (l,m) and (n,k) 
connected to the nodes to and n, respectively. Thus, the results 
of [13] and [6] no longer apply in the unbalanced setup. 
An analytical characterization of the flow region T in the 
unbalanced case is challenging because of the number of volt- 
age angles involved and the aforementioned coupling of the 
line power flows. Nevertheless, the following simple examples 
illustrate why one should expect a rank 1 solution from the 
relaxed OPF even in an unbalanced setup. 
Consider a 2-node unbalanced network, and suppose that a 
two-phase line connects the two nodes. Let Z 12 = [(0.0753 + 
jO.1181), (0.0156 + j0.0502); (0.0156 + j0.0502), (0.0744 + 
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jO.1211)] (as in [30]) and Y„ = 0, and assume that |V$| = 
1 for n = 1,2 and <fi = a,b. Fig. Ufa) depicts the feasible 
region 7"i 2 of powers (P^, iV+i) for 0? = 0°, 6»f 2 € 
[-180°, 180°], 0f G [110°, 130°], and 0f G [110°, 130°]. 
Notice that J 7 ^ is given by the Minkowski sum of the per- 
phase regions {{Pf^, P^i)}- It can be seen that the feasible 
flow region J-i 2 (the dark gray area) is a perturbed ellipsoid; 
specifically, expanding the expression of P m ^„, it follows 
that 0" b and 2 6 entail a perturbation of the center and of 
the axes of the ellipsoid that would be obtained if the line 
was balanced. Next, consider minimizing a strictly increasing 
function C(Pi^ 2 , fi->2) over .Fi 2 . It follows that the set of 
Pareto optimal points are the ones represented by the red dots 
in Fig. [TJa). Consider now the convex hull of T\2, which 
amounts to augmenting the dark gray area with the light gray 
one shown in Fig. |TJb). Clearly, points belonging to the light 
gray area lead to a solution of (P3) with rank higher than 1. 
Notably, the Pareto points of .Fi 2 and the ones of the relaxed 
region in Fig. |TJb) coincide; since C( ) is strictly increasing, 
the solution of (P3) must be on the Pareto boundary; thus, 
the optimal solution of (P3) has rank 1, and it is an optimal 
solution also for the nonconvex problem (P2). Granted that the 
Pareto points in Fig.sQJa) and (b) are the same, [13, Lemma 5] 
can be used to show that the two Pareto regions coincide also 
when constraints on the voltages are involved. 

The region of injected powers for a 3 -node network is exam- 
ined next. Let Z 12 = Z 23 = [(0.0753 + j'0.1181), (0.0156 + 
jO.0502); (0.0156 + j0.0502), (0.0744 + jO.1211)], and as- 
sume that \V$\ = 1 for all n = 1,2,3 and = a, b. 
Finally, let P„ := P£ + P%, with P$ the power injected 
at node n — 1,2,3 and phase <\> = a, b. The gray area 
depicted in Fig. |2fa) corresponds to the region of feasible 
power injections I for — tan _1 (3fi{[Z]~ 2 1 ]i j i}/3{[Z]~ 2 1 ]i j i}) < 
9% rn < tan- 1 ^^ 1 ]!,!}/^^ 1 ]!,!}) [13], and for the 
line-line angles 9f,9f,9f G [110°, 130°]; specifically, 
tan-^^IZ^ 1 ]!,!}/^^^ 1 ]^!}) « 58°. The angle differ- 
ence 9^ b is confined in the set [110°, 130°], as higher values 
are not likely to happen in practice [30]. Clearly, the gray 
region is nonconvex. If one considers minimizing a strictly 
increasing function of the injected powers {Pi,p2,P3} over 
1, it follows that the Pareto optimal points are the ones color 
coded red in Fig. |2a). It can be noticed that the Pareto region 
does not change if one takes the convex hull of 1; therefore, 
the solution of (P3) has rank 1, and it is an optimal solution 
for (P2). Next, fix the power injected (or absorbed) at node 
3, and consider minimizing a strictly increasing function of 
Pi and P 2 . The two-dimensional region of feasible powers Pi 
and P 2 is depicted in gray in Fig |2jb). Again, the Pareto front 
of the gray region and the one of its convex hull coincide. 

These examples suggest that the nonconvex OPF problem 
(P2) and its relaxed counterpart (P3) share their optimal so- 
lution when -t&n-\K{[Z m \]^}/^{[Z-^ n ]^}) < 9 a nm < 
tan ^{[Z-^^/S^Z"* ]^}), and the angle differences 
between conductors Offi are small enough. This further mo- 
tivates efforts toward analytical characterization of the power 
injection region in unbalanced distribution systems. 



IV. Distributed Solution 

Albeit polynomial, the computational complexity incurred 
by standard interior-point solvers for SDP does not scale well 
with the number of nodes N, and the number of constraints 
N c [17]. This lack of scalability is incurred also by alter- 
native methods based on off-the-shelf solvers for nonlinear 
programs [10]. Furthermore, the communication overhead 
required to collect data from all end users at the MEM, 
and subsequently disseminate the OPF solution back to the 
LCs, may lead to traffic congestions and substantial delays 
in the data delivery. Therefore, a distributed SDP solver, 
with minimal computational and communication costs, is well 
motivated. A distributed approach is also desirable in order to 
address possible concerns regarding data privacy and integrity, 
and when the microgrid includes single- or multi-facilities that 
are managed independently from the rest of the network in 
order to pursue specific economic interests [18]. 

Consider partitioning the microgrid into L areas {.4^ C 
J\f}fLl, and suppose that each area is controlled by an LC. 
In the distributed SDP solver to be derived, each LC will 
solve an OPF problem of reduced dimension for its controlled 
area. Let AS 1 ' be an "extended" area defined as A^ := 
A {l) U {n\(m,n) G £,m G A {l \n e A {i) ,1 ^ j}; that 
is, collects also the nodes belonging to different areas 
that are connected to A^ by a distribution line. Sets {^O} 
can be interpreted as counterparts of the "regions" considered 
in [20] in the context of transmission systems. Based on 
the overlaps among {A {l) }f =1 , define the set of neighboring 
areas for the Z-th one as MT® := {j\A {l) n A u) ^ 0}. 
Finally, let the vector v' z ' stack the complex line-to-ground 
voltages of the nodes in A^ 1 ' (that is, {v„}„ e _4 ( ()), and let 
*&,*m->n *nd V« := denote 
the sub-matrices of <frp„, 3?q n , *y„, 3>m^n and V, respec- 
tively, formed by extracting rows and columns corresponding 
to nodes in A^ 1 '. With these notational conventions, it is 
possible to re-write the SDP (P3) as 

L 

(P4) min 5^C$?(VW) (15a) 

2 = 1 

s.t. V w G B (l) , l = l,...L (15b) 

v y o (15c) 

where denotes the set of sub-matrices satisfying the 
following constraints per area A"': 

Tr(*£;' n V (, >)+P£„ = 0, V </>, Vn e A {1) nS (16a) 
Tr^VW) + <„ - ^,„Tr(*^VW) = 0, 

V0, Vne A {l) \(A {l) r\S)(16b) 
pgta < Tr(#^V ( ')) + P* >3 < P£ ,a *, V 0, Vn G «S« (16c) 
< Tr(*^VW) + Ql s < V 0, Vn G S (16d) 
(C°) 2 < Tt(*ty n V®) < (Kr x ) 2 ,V 0, Vn G .4 (() (16e) 

with the additional constraint {V"']v ,Vo = v o v o^ f° r tne area 
that contains the PCC. Further, Cm(V^) represents the cost 
associated with the area A"'; that is, it collects the terms 
in (|9at (if m = 1) or (|9at (if m = 2) that pertain to A^ . For 
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example, if one wishes to minimize the distribution losses, the 
corresponding cost becomes 

<7f(V«))= £ Tr(*UnV«) 

(m,n)e£|m,ne^tW 

+ J] ^Tr(*^„V«) (17) 

with *'„ OIl := *' m ^„ + *L_^ TO - On the other hand, the 
expression for C\ •'(V*-^) is not unique, depending on the 
specific network topology and operational setup. For instance, 
if ones assumes that each area is formed by one lateral or 
one sub-lateral [18], then C^(V^) accounts for the cost of 
power flowing into the (sub-)lateral from the "backbone" of 
the microgrid, and the cost of power generated within the 
(sub-)lateral. Instead, if customers are allowed to use DG 
units only to satisfy their own needs, C$\vW) boils down 

to C»(VW) = c„Cf(V«) + L^.n S c s E^4. 
Either way, the equivalent formulation (fT8l effectively ex- 
presses the cost as the superposition of local costs, and divides 
network constraints on a per-area basis. However, even with 
such a decomposition the main challenge lies in the PSD 
constraint (I15cb that couples local matrices {V^}. Indeed, 
if all submatrices {V^} were non-overlapping, the PSD 
constraint on V would simplify to y per area I, and 
(P4) would be decomposable in L sub-problems. However, 
this equivalence fails to hold here, since submatrices share 
entries of V. The idea is to identify valid network topologies 
(that is, valid partitions of the microgrid in smaller areas) for 
which the PSD constraint decomposition is feasible. 
To this end, define the following two auxiliary graphs: 

i) a "macro" graph Qa, where nodes represent the areas and 
edges are defined by the neighborhood sets {A/ - ^}^; and, 

ii) a "micro" graph Qm induced by the sub-matrices 

{y }l=i> that i s ' a S ra ph w i m 12 n eN l^™! n °d es ( one P er 
phase and node), with an edge connecting the nodes represent- 
ing the voltages V$ and Vj^ if the entry of V corresponding to 
V*(y*)* is contained in one of the sub-matrices {V"'}^. 
Examples of macro graphs are provided in Figs. |3]and|4] 

Based on these auxiliary graphs, results on completing 
partial Hermitian matrices will be leveraged to obtain PSD 
ones [24]. These results rely on the so-termed chordal property 
of the graph induced by {V' ! '}^ =1 to establish the equiv- 
alence between positive semidefiniteness of the overall matrix 
V and that of all submatrices corresponding to the graph's 
maximal cliques. Towards decomposing the PSD constraint 
into local ones, the following assumptions are made, which 
naturally suggest valid partitions of the microgrid: 
(Asl) The graph Qa is a tree; and, 

(As2) \A®\(A® r\A®)\ > for all i,l = l,...,L, i ^ I; 

that is, no nested extended areas are present. 
Condition (Asl) is quite reasonable in tree distribution net- 
works; for example, an area can be formed by a pair of nodes 
that are connected by a distribution line [6], or by laterals 
and sub-laterals [18]. (As2) is a technical condition ensuring 
that the subgraph induced by V' ! ' is a maximal clique of 
Qj\f, which allows using the results of [24]. Based on these 
assumptions, the following can be readily proved. 



Proposition 2: Under (Asl) and (As2), the graph is 
chordal, meaning that each of its cycles comprising four or 
more nodes has a chord. Furthermore, all its maximal cliques 
correspond to the elements of {V«}L. □ 

As established in [24], the PSD matrix V is "completable" 
if and only if Qjy is chordal, and all its submatrices corre- 
sponding to the maximal cliques of Qj^ are PSD. Therefore, 
constraining V to be PSD is tantamount to enforcing the 
constraint on all local matrices >: 0, Wl = 1, 

Notice that (Asl) requires the macro graph Qa to be a tree, 
while no conditions are imposed on the topology of the 
microgrid. Thus, it may be possible to find network partitions 
with an associated chordal graph also in the case of 
weakly-meshed microgrids. 

Next, let Vij := U»e(.#0n,Stf)) fn collect the indexes 
corresponding to the voltages {{V^}^ & -p n } of the nodes that 
the extended areas A^ 1 ' and A^ share. For example, if areas 
1 and 2 share nodes n = 5 and n = 6, then V12 indexes the 
voltages {Vjf},jJe7>5 an d {Ve}<pev 6 - Further, define as 
the submatrix of V") collecting the rows and columns of V^) 
corresponding to the voltages in Vij. With these definitions, 
and assuming that (Asl) and (As2) hold, problem (P4) can 



be re-written in the following equivalent form: 

L 

(P5) min VcW(VW) (18a) 
{v(0} ^ 

s.t.V (,) e6 (l) , l = l,...L (18b) 

vf=Vp\ l = l,...,L (18c) 

V (i) ^0, l = l,..,L (18d) 



where constraint (118ct enforces neighboring areas to consent 
on the entries of V") and V"' that they have in common. 
Clearly, constraints (118c) couple the optimization problems 
across areas. To enable a fully distributed solution, con- 
sider introducing the auxiliary variables {W^ IjeArco and 
{X^- }j 6 jv-(o per area. With these auxiliary variables, (P5) 
can be equivalently re-stated as 



L 

>) min y 

{V('> t O}^ 


C«(VW) 






(19a) 


s.t. V» 6 B« , 


1 = 1,. ..L 






(19b) 




jetf®, l = 


1,- 


..,L 


(19c) 


3{vf} = xp } , 


j e AT®, 1 = 


1, . 


..,L 


(19d) 


wp } = wf , 


ieJfC, l = 


l,- 


..,L 


(19e) 


x«=x», 


jetf®, 1 = 


l, ■ 


..,L. 


(19f) 



A similar approach was followed by [19], which utilized 
either primal or dual iterations to distribute the OPF in bal- 
anced transmission networks. A distributed OPF for balanced 
distribution feeders was derived in [6], where the dual (sub- 
gradient ascent was used. Unfortunately, sub-gradient ascent 
methods do not always lead to a satisfactory solution; when the 
dual function is non-differentiable and the step size is fixed, 
dual and primal iterates converge only on the average. What 
is more, recovering the primal variables from the optimal dual 
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variables is not always guaranteed [31, Sec. 5.5.5]. Besides 
dealing with unbalanced power networks, the novelty here 
consists in solving (P6) distributedly by resorting to the 
ADMM [25, Sec. 3.4], a powerful scheme that has been 
successfully applied to distributed optimization and estimation 



[SI'] Update V«(« 
( P7 (0) 



1) per area I = 1, . . . , L as: 



in several contexts. To this end, let {if} and {A-^ 1 '} be 
the multipliers associated with dl9cfc and dl9dt , respectively, 
and consider the partial quadratically-augmented Lagrangian 
of (O as 

L 

£({V (i) }, {Wf }, {Xf }, {if}, {A.f }) = £ {^'(V (i) ) 
+Tr((A< 



.CO 



min £C0( V (0 ){a £ }) 
vWho,{«i>o,ft>o} 


(25a) 


s.t. V» e B« 




(25b) 




^ I 


r< 0, Vj g A/"' 7 - 1 


(25c) 




"ft 


d o, Vj g M {1) 


(25d) 



where the local Lagrangian (125 at is given by 



if ) r (H{Vn -W[")) + SPOT} - W 



+ 

jeA/C) 



r 0h 



£ (l) (v (i) ,{^,ft}) 



3{V"}-X[ 



(j) 



+ ps*OT} 



x«||| 



]} (20) 



where k G K + is a positive constant [25, Sec. 3.4]. Then, the 
ADMM amounts to iteratively performing the following steps 
(i denotes the iteration index): 

[SI] Update primal variables: 

{V«(< + 1)}:= 

arg min £({V«}, {W« (i)}, {Xf (i)}, {if (i)}, {A« (i)}) 



+ £ [Tr((rf(i))^{Vf})+Tr((Af(i))^{V»})] 

and the vectors r 3 sr and r^g collect the real and imagi- 
nary parts, respectively, of the entries of the matrix V 



(0 



V 



(J), 



, « + vp } (i) 

[S2'] Update dual variables locally per area 1 = 1, 



-.(0 



(i + 1) = if (i) + £ (»{v« (< + 1)} - R{v«> (z + 1)) 



(26) 



s.t. V (0 G B w , 2 = 1 



CO 



.L. 



(21) Af(i + 1) 



(<) + 



{Vf (i + l)}-9f{V^(i + l) 



(j). 



[S2] Update auxiliary variables: 



{wf(i+l),Xf(i + l)} 



(0/ 



arg 



min £({V l ' ; (i + 



{w (0 x (0 

L 3 1 3 



} 



(0 (* + 1)}, {w«}, {x« } , {if (*)}, {AfMlivWCi), 



Furthermore, for any 



> 



(27) 



the iterates 



lf } (i),Ay ; (i)} produced by [S1']-[S2'] are 
convergent, and lim^ +00 V"C0 (i) = vf t for all I = 1, . . . , L, 



s.t. W ; u; = W 



(0 



(22) 



of 



[S3] Update dual variables: 



if (i + 1) = If (i) + K (K{V w (i + 1)} - W\ ] >(i + 1)) (23) 



to, 



with {V^pj} sub-matrices of the optimal solution V opt 
(P3). 

Proof. See the Appendix. □ 



Af (t + 1) 



Af (i) + K (3{V w (i + 1)} - X ; u; (i + 1)) (24) 

In step [SI], the per-area matrices {"V"'(i)} are obtained by 
minimizing (EQj, where variables {wf\i + l),xf\i + 1)} 
and the multipliers are kept fixed to their previous iteration 
values. Likewise, the auxiliary variables are updated in [S2] 
by fixing {V^(i + 1)} to their up-to-date values. Finally, the 
dual variables are updated in [S3] via dual sub-gradient ascent. 

Interestingly, the ADMM iterations [S1]-[S3] can be sim- 
plified by exploiting the favorable decomposability of the 
Lagrangian. To this end, the following lemma is first needed. 

Lemma 3: If the multipliers are initialized as if (0) = 



U), 



A«(0) 







\-P t3 \ 



then for every pair of neighboring areas 



I and j it holds that if (i) + T^' (i) = (W.I for each 



i > 1. Likewise, A^(i) + A\ J '(i) 



Af(0)=A^(0) = 0| 7 , i .|. 
Proof. See the Appendix. □ 

Using Lemma [3] steps [S1]-[S3] can be simplified as 
follows. Furthermore, convergence to the solution of the 
centralized SDP (P3) is established. 

Proposition 3: If if (0) = A^(0) = Q\ Vl .\, then [Sl]- 
[S3] boil down to the following primal-dual updates: 



.(0/ 



■Pi A 

0|7>,.| for each i > 1 if 



At each iteration, the LC of area I receives from the LCs 
of its neighboring areas j G A"t0 matrices Vj (i), and 
updates the local multipliers {if (i), A^(i)} via (I26b-(|27T>. 
These multipliers are locally stored at area I, and they are not 
exchanged among LCs (in contrast, multipliers are exchanged 
per iteration in [19]). Then, LC I updates VCO by solving 
(P7C0), and transmits Vj- (i + 1) to its neighboring areas 
j G A/"C0 . 

V. Numerical Tests 

The SDP-based solver is tested here using the following two 
networks, operating in a grid-connected mode: 

• the IEEE 37-node test feeder [30] shown in Fig. [3j and 

• the 10-node 3 -phase network depicted in Fig. [4] 

The optimization package CV>H, along with the interior- 
point based solver SeDuMiQ are employed to implement the 
centralized and distributed solvers in MATLAB. 

The 4.8 kV network of Fig. [3] is an actual portion of 
power a distribution network located in California; all the 
demanded complex powers are "spot" loads, and the network 
loading is very unbalanced [30]. Compared to the original 

3 [Online] Available: http://cvxr.com/cvx/ 

4 [Online] Available: http://sedumi.ie.lehigh.edu/ 
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Fig. 3. First Test: IEEE 37-node feeder. 
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Fig. 5. Power generated by the DG units [kW]. 

The minimum and maximum utilization and service voltages 
are set to V r f n = 0.95 pu and V„ max = 1.05 pu for all nodes. 
Thus, voltage regulation is enforced without requiring changes 
in the voltage regulator taps (as in e.g., [15]). Further, the 
voltage at the PCC is set to v a = [1Z0°, 1Z-120°, 1Z120°] T 
pu. The average computational time required by SeDuMi to 
solve the centralized problem (P3) was 9.0 sec and 0.3 sec 
(machine with Intel Core i7-2600 CPU @ 3.40GHz), which 
is significantly lower than the time required by commercial 
solvers for non-linear programs (see e.g., [10]). 

TABLE II 

Results for the IEEE 37-node feeder 



Cs/C 


floss [kW] 


C*2 [$] 


P [MW] 


P G [MW] rank(V opl ) 





39.30 


57.9 


1.4463 


1.0500 1 


0.25 


39.30 


68.4 


1.4463 


1.0500 1 


0.50 


39.30 


78.9 


1.4463 


1.0500 1 


0.75 


38.97 


89.3 


1.4469 


1.0491 1 


I 


36.60 


99.7 


1.5979 


0.8957 1 


1.25 


57.22 


102.6 


2.3088 


0.2054 1 



Fig. 4. Second test: (a) 10-node microgrid; (b) corresponding graph CI- 



TABLE III 

Results for the 10-node network in Fig. [4] 



scheme however, 7 DG units are placed at nodes S = 
{10,12,16,19,26,32,36}. Specifically, single-phase conven- 
tional DG units supply a maximum real power of 50 kW per 
phase, and are operated at a unit power factor (PF); that is, 
Qq" s = <9™s = for all seS. Line impedances and shunt 
admittances are computed based on the dataset in [30]. Finally, 
delta-wye conversions are performed whenever necessary. As 
for the network of Fig. Ha), the line admittances are all set 
to Z m „ = [0.0693 + jO.2036, 0.0312 + j0.1003, 0.0316 + 
j0.0847; 0.0312 + jO. 1003, 0.0675 + j0.2096, 0.0307 + 
j'0.0770; 0.0316 + j0.0847, 0.0307 + j0.0770, 0.0683 + 
j0.2070] n (see [30]), which gives rise to an unbalanced 
operation of the network. Shunt admittances are neglected. 
Single-phase conventional DG units are placed at nodes S = 
{5, 7}, and they can supply a maximum real power of 50 kW 
per phase, at unit PF. All the loads are "spot", and the loading 
is assumed balanced. Specifically, the active and reactive 
loads are set to {0,50,0,130,130,110,110,0,90,90,} kW 
and {0, 20, 0, 82, 82, 60, 60, 0, 30, 30, } kVAr, respectively, on 
each phase. The voltage magnitude at the PCC is 4.16 kV. With 
these two choices, the performance of the proposed approach 
can be assessed for different network sizes, line characteristics, 
loading, and different network partitions shown next. 



Cs/C0 


fl oss [kW] 


C 2 [$] 


P [MW] 


P G [MW] rank(Vopi) 





18.38 


59.9 


1.4984 


0.3000 1 


0.25 


18.38 


62.9 


1.4984 


0.3000 1 


0.50 


18.38 


65.9 


1.4984 


0.3000 1 


0.75 


18.38 


68.9 


1.4984 


0.3000 1 


1 


18.27 


71.9 


1.5457 


0.2526 1 


1.25 


23.08 


72.1 


1.8031 


0.0000 1 



Tables HI] and [TTT1 list the real power drawn at the PCC 
Pq := J2(f>£-p Vo(Io)*' tne total power generated by the DG 
units Pq := J2 s P^, s , and the overall power losses Pi oss 
and costs of supplied power, when the cost (O is employed. 
The costs of supplied power are set to cq = 40 $/MW and 
ci = ... = c 7 e {0,10,20,30,40,50} $/MW. Notice that 
minimizing (O is tantamount to minimizing the power loss (Q]) 
when Co = 1 and c s = 1 for all s € S. In fact, it can be clearly 
seen in Tables HI] and HII] that the power loss is minimized for 
this choice of {c s }. Powers Pi oss , Po, and Pa remain the same 
when c s < cq\ however, when c s > cq (which holds during 
the late-night hours) the DG units reduce the generated active 
powers; further, the power loss becomes significantly higher 
when DG units are not used at the maximum extent. Fig. |5] 
depicts the active powers generated by DG units as a function 
of c s /cq, for the IEEE 37-node test feeder. It can be noticed 
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Fig. 6. Convergence of the ADMM, with feeder partitioned as in Fig. [5] 
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Fig. 7. Convergence of the ADMM, network partitioned as shown in Fig. [4] 

that the DG units electrically close to the PCC are not utilized 
when c s > co; on the other hand, DG 6 and 7 still operate at 
more than 50% of their maximum capacity. 

Interestingly, the rank of matrix V opt was always 1 . There- 
fore, the globally optimal solutions of (P2) (and hence of 
the nonconvex (PI)) were always attained. In other words, 
no lower power losses or costs of supplied power can be 
attained with alternative OPF solution approaches. The rank 
of matrix V op t was greater than 1 for the IEEE 37-node test 
feeder when c s /cq > 1.75 (two non-zero eigenvalues, with 
Ai oc 10 2 and A2 oc 10~ 2 ). Nevertheless, a rank-1 solution 
was readily obtained upon raising the voltage magnitude at 
the PCC from 1 to 1.02 pu. This further prompts an analytical 
characterization of the feasibility region of (PI). Additional 
tests were performed on the IEEE 13-node feeder [30], and 
rank-1 matrices V opt were again always obtained (results are 
not reported here due to space limitation). 



Convergence of the proposed distributed SDP solver is 
showcased in Fig. [6] where the 37-node feeder is partitioned 
as shown in Fig. [3] This partition resembles the case where 
laterals include multi-facilities that are managed independently 
from the rest of the network [18]. The trajectories correspond- 
ing to AV Lj (i) := ||V; (j) (i) - V^(i)||i/36 per iteration i, 
are reported for different values of the ADMM parameter k, 
and are compared with the ones obtained by using the sub- 
gradient ascent-based distributed algorithm developed in [6] (s 
denotes the step size, which is assumed to be either constant, or 
monotonically decreasing). It can be noticed that the proposed 
distributed solver exibit a considerably faster convergence 
than the one based on the sub-gradient. Considering that 
the nominal voltage of the feeder is 4.8 kV, in less then 
than 50 iterations the average gap between the entries of 
Vp } ( Z ) and Vf\i) is on the order of a few volts. Further, 
notice that the convergence rate is approximately the same for 
AVi.2(i), AVi.3(i), and AVi^(i). The average computational 
time required by SeDuMi to solve each sub-problem was 
of 6,3,0.3, and 2.7 sec. Lower computational times can be 
obtained by selecting areas of smaller size. 

Fig. [7] illustrates the convergence of the distributed SDP 
solved when the considered 10-node network is partitioned as 
shown in Fig |4] The ADMM-based method outperforms the 
one based on the sub-gradient here too. In this case, the gaps 
between the elements of matrices {V> (i)} rapidly vanish 
when k = 100 and k = 10 after approximately 100 iterations. 
In this case, the average time required by SeDuMi to solve 
each sub-problem was 0.1 sec. 

VI. Concluding remarks 

The OPF problem was considered for microgrids oper- 
ating in an unbalanced setup. Inspite of the inherent non- 
convexity, the SDP relaxation technique was advocated to 
obtain a convex problem. As corroborated by numerical tests, 
the main contribution of the proposed approach consists in 
offering the potential to obtain the globally optimal solution 
of the original nonconvex OPF. A distributed SDP solver was 
also developed by resorting to the ADMM. The distributed 
algorithm ensures scalability with respect to the microgrid 
size, robustness to communication outages, and preserves data 
privacy and integrity. 

Appendix 

Proof of Lemma\\\ Equality (TTab can be readily established 
by noticing that \V*\ 2 = v w e£(e^) r v = Tr(e^(e*) r V). 
To prove (fTbl i. notice first that the injected apparent power at 
node n and phase (f> is given by V*(I$)* = (V$*I%)* = 
(y H ef l {ef l ) r \) n . Next, using i = Yv, it follows that 
(v*e£(e£)Ti)K = (v«ef(ef)^Yv)« = (v*Y*v)* = 
v w (Y^) w v, which can be equivalently rewritten as Tr(Y^V). 
Thus, the injected real and reactive powers can be obtained 
by using, respectively, the real and imaginary parts of (Y*) w . 
Finally, to prove d7dl i, define first the |P mra | x 1 vector i m ~m 
collecting the complex currents flowing from m to n on each 
phase, and notice that i m ->n = Z~^(v m - v„) = A,„^„v 
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and v m = B m v. Then, it follows that P m _>„ = 5R{i„^„v} = 
5R{v«A^„B m v} = K{Tr(A^„B m V)}. 

Proof of Lemma [3] The proof is provided for the multipliers 
r^. 1 , same steps can be followed for the dual variables A^. . 

Suppose that i = 1, and notice that the optimiza- 
tion problem ( 1221 decouples per pair of neighboring ar- 
eas; specifically, a number of sub-problems is to be 
solved with respect to (wrt) the only pair of variables 
(W ; (3) ,wj ) for j e AfW (and clearly I € AfW). Ne- 
glecting irrelevant terms, and setting the first-order deriva- 
tive of £({VW(l)},{wf },{xf},{lf (0)},{Af (0)}) 
wrt (Wp^ , W^) to zero, it readily follows that the minimizer 
Wj- (1) is given by 

wf(l) = ^(vf (lj+vf^l)). ( 28 ) 

Next, substituting d2HJ into [S3], and setting 1^(0) = 
rf(0) = yields rf(l) = § (Vf (1)-V, W) (1)) and 
lf°(l) = § (Vj 0) (l) - vj°(l)Y It then follows by induc- 
tion that r^(i) = — rp^(i) also for all subsequent iterations. 

Proof of Proposition [3] Substituting (|28l > into [S3], and 
setting 1^ (0) = I 1 ^ (0) = 0, one can readily obtain (l26t- 
(|27] |. Consider now [SI']- First, substitute d28j into d20t . and 
discard irrelevant terms. Next, introduce an auxiliary variable 
ctj to upper bound the squared I2 norm of rj.sft, and add 
the constraint 1 1 r ^ , gpj 1 1 2 — a j- Thus, (I25dl ) follows by Schur's 
complement. Indeed, it is certainly possible to introduce one 
single variable a to upper bound the squared I2 norm of 
the vector [r^, r^ 3 , . . . , r^ (i)kjR , rJ m) ^] T , and derive a 
linear matrix inequality similar to (125 dt by using Schur's 
complement. Finally, since the cost in (P6) is convex, and 
the constraint set is bounded and convex, every limit point 
of {V'''(i)} is an optimal solution to (-P6), as established 
in [25, Prop. 4.2]. Since (P6) and (-P3) are equivalent under 
(Asl) and (As2) hold, every limit point of {V' ! '(i)} is also 
an optimal solution to (P3). 
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